Revisiting the effect of external fields in Axelrod's model of social dynamics 
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The study of the effects of spatially uniform fields on the steady-state properties of Axelrod's 
model has yielded plenty of controversial results. Here we re-examine the impact of this type of field 
for a selection of parameters such that the field-free steady state of the model is heterogeneous or 
multicultural. Analyses of both one and two-dimensional versions of Axelrod's model indicate that, 
contrary to previous claims in the literature, the steady state remains heterogeneous regardless of 
the value of the field strength. Turning on the field leads to a discontinuous decrease on the number 
of cultural domains, which we argue is due to the instability of zero-field heterogeneous absorbing 
configurations. We find, however, that spatially nonuniform fields that implement a consensus rule 
among the neighborhood of the agents enforces homogenization. Although the overall effects of the 
fields are essentially the same irrespective of the dimensionality of the model, we argue that the 
dimensionality has a significant impact on the stability of the field-free homogeneous steady state. 

PACS numbers: 87.23.Ge, 89.75.Fb, 05.50.+q 



I. INTRODUCTION 

Axelrod's model of social dynamics was introduced to 
explore the mechanisms behind the persistence of cul- 
tural differences in a society The agents are repre- 
sented by strings of cultural features of length F, where 
each feature can adopt a certain number q of distinct 
traits. The interaction between any two agents takes 
place with probability proportional to their cultural sim- 
ilarity, i.e., proportional to the number of traits they have 
in common. The analysis of this model by the statisti- 
cal physics community has revealed a rich dynamic be- 
havior with a nonequilibrium phase transition separating 
the culturally heterogeneous from the culturally homoge- 
neous regime [2|-|4|. 

An interesting characteristic of Axelrod's model, which 
sets it apart from most lattice models that exhibit 
nonequilibrium phase transitions @, is that all station- 
ary states of the dynamics are absorbing states, i.e., the 
dynamics always freezes in one of these states Q. In fact, 
according to the rules of Axelrod's model, two neighbor- 
ing agents who do not have any cultural trait in com- 
mon cannot interact and the interaction between agents 
who share all their cultural traits does not result in any 
change. Hence at the stationary state we can guarantee 
that any pair of neighbors are cither identical or com- 
pletely different regarding their cultural features. This 
allows us to easily identify the stationary regime, which 
is a major problem in the characterization of nonequi- 
librium phase transitions The problem, however, is 
that the dynamics can take a very large time to freeze to 
a homogeneous configuration for some initial conditions 
@> 111 ) which is the main reason there are so few numeri- 
cal estimates of the transition lines of the phase diagram 
of Axelrod's model [4]. 

From the perspective of the statistical physics, the ap- 
pealing feature of Axelrod's model in a lattice of dimen- 
sion d is the existence of a threshold value q c = q c (F) 
with F > 1 below which the stationary regime is mono- 



cultural (i.e., spatially homogeneous) and above which 
multicultural (i.e., spatially heterogeneous). This result 
holds true for both the two-dimensional @, 0] and the 
one-dimensional [7[ versions of Axelrod's model. We re- 
call that the sources of disorder in this model are the 
stochastic update sequence and the choice of the initial 
configuration: it is the competition between the disorder 
of the initial configuration and the ordering bias of the 
local interactions that is responsible for the nontrivial 
threshold phenomenon. 

The introduction of an external global field to influ- 
ence the agents' beliefs aiming at modeling the effect of 
the mass media Q resulted in a surprisingly difficult and 
controversial problem since the external field favored the 
heterogeneous instead of the homogeneous regime as one 
would naively expect (see Q for an explanation of this 
finding). In addition, a considerable amount of effort has 
been devoted to searching for a threshold on the intensity 
of the media influence such that above that threshold, the 
community is multicultural and below it, the community 
is mono-cultural (see, e.g., We have shown, 

however, that this threshold is an artifact of finite lat- 
tices, and that even a vanishingly small media influence 
is sufficient to produce cultural diversity in a region of 
the parameters space where the homogeneous regime is 
dominant in the absence of the media [3, [HJ . 

In this paper we address another curious finding re- 
garding the effect of the media (or external field) in Ax- 
elrod's model: if the control parameters are such that 
the stationary regime is multicultural at zero field, then 
turning the field on will lead to a homogeneous state in 
the limit of vanishingly small field intensity (TlJ. Here 
we argue that this claim is not correct and that the mul- 
ticultural regime remains multicultural, though with a 
reduced cultural diversity, regardless of the intensity of 
the global external field. By global field we mean a field 
that is spatially uniform (i.e., it is the same for all agents) 
although not necessarily time independent @. 

Over and above the re-examination of far-out claims 
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on the effect of external fields in Axelrod's model, in this 
contribution we show that the effect of the media in the 
one-dimensional model is qualitatively identical to the 
two-dimensional model, which has been the sole focus of 
systematic analyses up to now (see Q for the study of 
the media- free one-dimensional Axelrod's model). Since 
simulations of the one-dimensional model are fast we will 
use their results as clues to expose the properties of the 
stationary state of the two-dimensional model in the com- 
putationally prohibitive regime of large lattices. 

The paper is organized as follows. In Section |TT] we 
briefly present the basic elements of Axelrod's model and 
describe the different types of external global fields (me- 
dia) we study in this paper. In Section lHll we show the re- 
sults of the simulations for the one-dimensional model in 
the two cases of interest, namely, when the media-free ab- 
sorbing configurations are homogeneous and when they 
are heterogeneous. In the first case, we show that, simi- 
larly to the two-dimensional version [j| [l4| , even a van- 
ishingly small field intensity is sufficient to break up the 
homogeneous steady state. In the second case we show 
that the external field reduces the number of cultural 
domains but the steady state remains heterogeneous, re- 
gardless of the field strength. This conclusion is corrob- 
orated by the simulations of the two-dimensional model 
described in Section HVl A different type of field - a spa- 
tially non-uniform field that implements the consensus 
or majority rule among the neighborhood of the agents 
- is discussed in a separate section because it cannot be 
interpreted as a media and produces results completely 
distinct from the spatially uniform fields (see Section fV|) . 
Finally, in Section IVII we summarize our main findings 
and present our concluding remarks. 



II. MODEL 

In the original formulation of Axelrod's model [![, 
which we will adhere to here, each agent is character- 
ized by a set of F cultural features which can take on 
q distinct values. In the two-dimensional version, the 
agents are fixed in the sites of a square lattice of linear 
size L with free boundary conditions (i.e., agents in the 
corners of the lattice interact with two neighbors, agents 
in the sides with three, and agents in the bulk with four 
nearest neighbors), whereas in the one-dimensional vari- 
ant the agents are fixed in the sites of a chain of length L 
with the same boundary conditions. From our perspec- 
tive, the advantage of using free boundary conditions is 
the easy of implementing Hoshen and Kopelman algo- 
rithm for counting the number of clusters in a lattice 
15], but since we are interested in the properties of the 
steady state for very large lattice sizes, the choice of the 
boundary conditions is largely irrelevant for Axelrod's 
model. 

The initial configuration is completely random with the 
features of each agent given by random integers drawn 
uniformly between 1 and q. At each time we pick an 



agent at random - the target agent - as well as one of 
its neighbors. These two agents interact with probabil- 
ity equal to their cultural similarity, defined as the frac- 
tion of common cultural features. An interaction consists 
of selecting at random one of the distinct features, and 
making the selected feature of the target agent equal to 
its neighbor's corresponding trait. This procedure is re- 
peated until the system is frozen into an absorbing con- 
figuration. Thus Axelrod's model can be viewed as F 
coupled voter models [l6|. 

The introduction of an external field or global media 
in the standard model follows the ingenious suggestion 
of adding a virtual agent which interacts with all agents 
in the lattice and whose cultural traits reflect the me- 
dia message [8(. In the original version, each cultural 
feature of the virtual agent has the trait which is the 
most common in the population - the consensus opinion. 
Henceforth we will refer to this type of external field as 
the consensus field. The second type of field we consider 
is constant in time, i.e., the media message is fixed from 
the outset, so it really models some alien influence im- 
pinging on the population. We will refer to this field as 
the static field. Explicitly, we generate the culture vector 
of the virtual agent at random and keep it fixed during 
the dynamics [lO, [H[ • (These two types of field were re- 
ferred to as global and external media by Ref. [ll|, but 
since both fields are global and external here we opt to 
use a more informative nomenclature.) 

Regardless of the type of external field, the interaction 
of the media (virtual agent) with the real agents is gov- 
erned by the control parameter p G [0, 1], which may be 
interpreted as a measure of the strength of the external 
field influence. As in the original Axelrod's model, we 
begin by choosing a target agent at random, but now it 
can interact with the media with probability p or with 
its neighbors with probability 1 — p. Since we have de- 
fined the media as a virtual agent, the interaction follows 
exactly the same rules as before. The media-free model 
is recovered by setting p = 0. 

Because of the unusually large times needed for some 
initial configurations to relax to a homogeneous absorb- 
ing configuration - relaxation to heterogeneous configu- 
rations is typically very fast - the simulation of the dy- 
namics of Axelrod's model must be made as efficient as 
possible. In our simulations we consider two lists of active 
bounds. The first list (list A) is composed by the active 
bounds that connect real agents, whereas the second list 
(list B) contains the active bounds that connect the vir- 
tual agent (media) with the real ones. In both cases, 
an active bound is defined as a bound that connects two 
agents that have at least one feature in common and at 
least one feature distinct from each other. Here instead 
of picking the target agent at random we first select one 
of the two lists - list A with probability 1 — p and list B 
with probability p - and then pick a bound at random 
from the selected list. In case of a bound from list B, the 
target agent is of course the real one but if the bound be- 
longs to list A we choose the target agent at random from 
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the two options. In the case that the cultural features of 
the target agent are modified by the interaction with its 
neighbor, we need to re-examine the active/inactive sta- 
tus of all bounds associated to the target agent so as to 
update the lists of active bounds. The dynamics is frozen 
when the two lists of active bounds are emptied. 



III. EXTERNAL FIELDS IN THE 
ONE-DIMENSIONAL MODEL 

In the attempt to make Axelrod's model more 'realis- 
tic', researchers have studied the model in a variety of 
complex networks (see, e.g., [I34i3)> w ^ n the usua l re- 
sult that the multicultural regime is destabilized by the 
increase of the connectivity of the network and so the 
stationary regime is homogeneous regardless of the val- 
ues of the parameters q and F. It is curious that the 
simple one-dimensional variant that preserves the phase 
transition Q received comparatively almost no attention. 
In this section we show that the effect of external fields is 
essentially the same in one and two dimensions as far as 
the favoring of the homogeneous or heterogeneous is con- 
cerned. The dimensionality introduces some distinctive 
effects, however, as we discuss next. 



A. Field-free homogeneous regime 

We begin our analysis with the once puzzling situation 
in which the homogeneous media-free stationary regime 
becomes heterogeneous under the influence of the exter- 
nal media Since this problem was extensively stud- 
ied in the two-dimensional lattice (9I-H1I H3 | we present 
here only a brief analysis of the effect of the static me- 
dia, aiming at highlighting the similarities and differences 
between the one and two-dimensional results. 

To characterize the steady state of Axelrod's model 
we focus on two basic statistical measures, namely, the 
normalized average number of domains Ndom/ L and the 
average relative size of the domains that do not belong 
to the largest domain 1 — S max /L. These quantities are 
shown in Fig.[T]togethcr with the fraction of runs trapped 
into homogeneous absorbing configuration £/, . A domain 
or cluster is a bounded region in which all agents share 
the same culture. We note that in this figure, as well as in 
the next figures of this paper, the statistical error bars are 
smaller or at most equal to the symbol sizes. Typically, 
each symbol represents the result of the average over 10 3 
to 10 4 independent runs of the stochastic dynamics. 

The first point to note is the remarkable similarity 
between the results presented in the upper panel of 
Fig. [T] for the one-dimensional and those for the two- 
dimensional model @ . For a fixed finite value of L there 
seems to exist a threshold value for the media strength 
p = p c below which the regime is mono-cultural (ToL fllj | . 
The quantity £,h shown in the middle panel illustrates 
this somewhat odd predominance of the homogeneous 
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FIG. 1: Results for the static media in the one-dimensional 
model with F = 5 and q — 3 and for media strengths (unfilled 
symbols, bottom to top in the upper and lower panels and 
top to bottom in the middle panel) p = 0.01,0.02,0.03,0.04 
and 0.05. The filled triangles are the media-free (p — 0) 
results. Here (Ndom) is the average number of domains, £j, is 
the fraction of runs that ended up in homogeneous absorbing 
configurations, and (S moJ; ) is the average size of the largest 
domain. The straight line in the upper panel is 1/L. 
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absorbing configurations, in which the agents are identi- 
cal to the static media, for lattices of intermediate size 
and small p (the case p = is discussed below). How- 
ever, as illustrated in the figure this 'threshold' decreases 
with increasing L and so it is a finite size effect. Thus 
our conclusion is that in the thermodynamic limit even a 
vanishingly small field is sufficient to break up the mono- 
cultural regime. More pointedly, extrapolating the ra- 
tio {Ndom)/L (upper panel of Fig. [T]) for L —> oo and 
plotting the result against p in a log-log graph yields 
limL^-oo (Ndom) / L ~ V 2 ' 1 in the limit p -> 0. 

The middle and lower panels of Fig. Q] reveal 
some remarkable differences between the one and two- 
dimensional models. Firstly, for p = only about 35% of 
the samples (random initial configurations) ended up into 
strictly homogeneous absorbing configurations, whereas 
for the two-dimensional model this happens for all sam- 
ples in the homogeneous regime. The reason we keep 
referring to this regime as the homogeneous regime is 
that the order parameters Nd om /L and S max /L take on 
values compatible with a homogeneous phase. This only 
happens because the heterogeneous absorbing configura- 
tions are composed of a single macroscopic domain to- 
gether with a non-extensive number of microscopic do- 
mains. Hence the probability that a randomly chosen site 
in such configuration belongs to the largest domain is 1 
in the thermodynamic limit. Secondly and more impor- 
tantly, in the two-dimensional model the fraction of the 
sites that are not part of the largest domain (the media 
in that case) is finite and approaches zero forp-j> @, 
whereas in the one-dimensional model that fraction tends 
to unity regardless of the field strength p > 0. This means 
that the effect of the field in one-dimension is quite ex- 
treme - a vanishingly small field is enough to completely 
destroy the uniform regime, producing a fragmented con- 
figuration composed of microscopic domains. (We know 
that the domains are microscopic because (S max )/L — > 
in the thermodynamic limit.) This result is probably a 
consequence of the 'weakness' of the uniform regime at 
p = discussed above. 
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FIG. 2: Normalized number of domains in the one- 
dimensional model for the static (upper panel) and consensus 
(lower panel) media as function of the chain size L. The filled 
symbols are for p = 1 (circles) and p = (triangles) whereas 
the unfilled symbols are for (top to bottom at L = 10 5 ) 
p = 0.99,0.6,0.4,0.2 and 0.01. The other parameters are 
F = 5 and q = 10. 



B. Field-free heterogeneous regime 

We turn now to the problem that motivated this pa- 
per, namely, the claim that the effect of a vanishingly 
small external field is to turn the field-free heterogeneous 
steady state into a field-induced homogeneous state in the 
limit L — > oo [ll[ . This result would then be analogous to 
the homogenizing effect of an external noise that changes 
traits at random with some small probability [2(| ■ Alas, 
we find no evidence of such mind-boggling effect neither 
in the one nor in the two dimensional Axelrod's model. 
Actually, the original small-lattice simulations of Ref. [ll[ 
do not seem to support the existence of such phenomenon 
either (see Sect. ITV|) . 

In Fig. [2] we show the average number of domains for 
both static (upper panel) and consensus (lower panel) 



media. We use filled symbols to display the results for 
the two extreme values of the field intensity (p = and 
p = 1), and unfilled symbols for the intermediary field 
strengths. This figure exhibits many noticeable results. 
The number of domains is maximum in the case only 
interactions with the field are permitted (p = 1). Al- 
lowance of local interactions between neighboring agents 
results in less heterogeneous configurations, as expected. 
The surprise is that, even for finite L, the average number 
of domains jumps to a lower value as p departs infinitesi- 
mally from unity (the same phenomenon happens in the 
two-dimensional model 11]). Although (Ndom) exhibits 
a smooth dependence on p S (0, 1) we observe another 
jump (now to a higher diversity value) at p = 0. The 
fact that these jumps occur at any finite value of L im- 
plies that the heterogeneous absorbing configurations are 
unstable to single-site changes. This conclusion is sup- 
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ported by the fact that the number of domains (and the 
diversity, as well) are reduced when the field or the local 
interactions are turned on. Put differently, the number of 
domains decreases (discontinuously) from p = to p > 
as well as from p = 1 to p < 1 . 

We note that for p = 1 the consensus media is fixed 
with the traits reflecting the consensus values for the ran- 
dom initial configuration of the agents. However, for fi- 
nite L the results differ somewhat significantly from those 
of the static media (see Fig. [2j. Why is that? The reason 
is that the traits of the static media are chosen randomly 
and independently of the also random initial traits of the 
agents, whereas in the case of the consensus media the 
media traits are not independent from the agents' traits, 
being given by the majority rule. This correlation has a 
strong effect for small chains but becomes irrelevant for 
large L since in this case the fraction of agents that share 
a given media trait is not much greater than l/q, which 
is the expected value of this fraction for the static media. 
In fact, the correlation between the media and the agents 
in the initial configuration explains the dips observed in 
the lower panel of Fig. [2] for p > since it decreases the 
number of agents with antagonistic traits with respect 
to the media, resulting in less fragmented configurations 
than in the case the media traits are set randomly. How- 
ever, as pointed before, this effect disappears for large 
chains. 

The discontinuities at p = and p = 1 are more easily 
visualized in Fig. [3] where we present the extrapolation 
to infinite chain sizes, g^ = lim l -> oo(-^V dom) / L, of the 
results exhibited in Fig. [2j The results for both types 
of media are practically indistinguishable. For p — > 
we find ffoo -> 0.38 ± 0.01 whereas g x -> 0.66 ± 0.01 
for p — > 1, which are valid for both media types (see 
the legend of Fig. [3] for the values of at p = and 
p = 1). The important point here is that the data offer no 
evidence whatsoever that g^ would vanish in the limit 
p — > as claimed by Ref. [ll| for the two-dimensional 
model. We should emphasize that the decrease on the 
number of domains induced by a vanishingly small field 
is the expected outcome of the assay, since the small field 
destabilizes some of the field-free domains, but lacks the 
strength to create new field-induced domains. 

In addition, we find that the average size of the largest 
domain (S ma x) grows as In L for large L (data not shown) 
whereas the size of a typical domain is on the order of 1, 
regardless of the media type and strength. These findings 
are similar to those reported for the majority- vote model 
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IV. EXTERNAL FIELDS IN THE 
TWO-DIMENSIONAL MODEL 

We turn now to the study of the two-dimensional Axel- 
rod's model, which is considerably more computationally 
demanding than the analysis of the one-dimensional ver- 
sion presented before. The effect of external fields on the 




FIG. 3: Extrapolation to L — > oo of the normalized number 
of domains in the one-dimensional model for the static (cir- 
cles) and consensus (triangle) media as function of the field 
strength p. At p = we find = 0.483 ± 0.001 and at 
p = 1 we find = 0.793 ± 0.001 for both media types. The 
parameters are F — 5 and q — 10. 



field-free homogeneous regime is well-understood by now 
[t| [3] : the results are similar to those exhibited in Fig. 
[T] except for the relative size of the largest domain which 
in the limit p — > tends to 1 in the two-dimensional case 
9] and to in the one -dimensional version (lower panel 
of Fig. [1]). Hence we will consider here only the effect 
of external fields on the field-free heterogeneous regime. 
Moreover, we focus most our efforts on the parameter set 
F = 2 and q = 8 rather than on the set F = 5 and q = 30 
of Ref. [ll|- In both cases the field-free (p — 0) regime 
is heterogeneous, but simulations using our parameter 
selection are much faster of course. 

Accordingly, in Fig. 2] we present the dependence of 
the average number of domains on the linear size of the 
square lattice. The results are qualitatively similar to 
those obtained for the one-dimensional model and sum- 
marized in Fig. [2] The correlation between the con- 
sensus media and the agents results in less heteroge- 
neous absorbing configurations in comparison with the 
configurations induced by the static media, but this dif- 
ference becomes negligible as the lattice size increases. 
This correlation effect is less dramatic than in the one- 
dimensional model because for the same value of L in 
the x-axis there are many more agents (L 2 to be pre- 
cise) in the two-dimensional lattice. The result of the 
extrapolation to infinite lattice sizes using the definition 
goo = lim l^oo(S dom) I L 2 are summarized in Fig. [SJ For 
p^Owe find Soo -s- 0.36±0.01 whereas g^ -> 0.53±0.01 
for p 1. As before, we observe discontinuities at p = 
and p = 1 which take place even for finite L, and no 
evidence at all that the g^ would tend to in the limit 
p — > as claimed by Ref. [ll|. We note, however, that 
extrapolation to L — > oo using lattices of linear size up to 
L = 70 as done in Ref. [Ill ] may result in a significant un- 
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FIG. 4: Normalized number of domains in the two- 
dimensional model for the static (upper panel) and consensus 
(lower panel) media as function of the linear size L of the 
square lattice. The filled symbols are for p = 1 (circles) and 
p = (triangles) whereas the unfilled symbols are for (top to 
bottom at L = 10 3 ) p = 0.99,0.4,0.2 and 0.01. The other 
parameters are F — 2 and q = 8. 



derestimate of g^ since, especially for small p, one may 
be misled by the transient region where gh = {Ndom)/L 2 
decreases abruptly with increasing L (see Fig. [4]). 

Our simulations for the selection of parameters F = 5 
and q = 30 used in Ref. [IH led to similar conclusions. 
As already pointed out the large relaxation times make 
an extensive analysis of this parameter set prohibitive, 
so we used a small number of samples (typically 100), 
which resulted in rather noisy data points. For small 
lattice sizes (L < 70) our results fully agree with those 
of [ll]. Actually, it should be said that those results do 
not support the claim that the vanishes in the limit 
p — > for the static and the consensus media (see Fig. 
7 of Ref. [H|). Perhaps this claim was motivated by an 
unwarranted generalization of the effect of a very distinct 
type of field as we discuss in the next section. 



FIG. 5: Extrapolation to L — > oo of the normalized number 
of domains in the two-dimensional model for the static (cir- 
cles) and consensus (triangle) media as function of the field 
strength p. The results are indistinguishable for the two me- 
dia types. At p = we find goa = 0.425 ± 0.001 and at p = 1 
we find <?oo = 0.869 ± 0.001. The parameters are F — 2 and 
9 = 8. 



V. LOCAL FIELD 

Here we consider a time and spatially nonuniform field 
introduced in Ref. [ll|, referred to as local field (or me- 
dia) . For a given target agent the traits of this field reflect 
the consensus trait of its nearest neighbors. So this field 
is a multi-state variant of the majority- vote rule [2ll - [2(| 
and, in that sense, we are reluctant to characterize it as 
a media, much less as a mass media. In fact, its local 
nature sets it apart from the static and consensus fields 
studied in the previous sections, which are spatially uni- 
form fields, and so we see no reason to assume the sys- 
tem will behave similarly under the effect of such diverse 
fields. 

The overly distinct effect of the local field on the steady 
state of Axelrod's model is illustrated in Fig. [6] for the 
one-dimensional variant. Most interestingly, in this case 
the steady-state is less heterogeneous when the agents 
interact with the field only (p = 1) than when the field is 
off (p = 0). This is a peculiarity of the one-dimensional 
model which suits well to illustrate the tendency to ho- 
mogenization of the local field. The results for all values 
of p G (0, 1) are indistinguishable from that exhibited in 
the figure for p — 0.99, which shows the dominance of 
homogeneous absorbing configurations already for small 
chains. In fact, for L > 50 the data points fall on the 
curve shown in the figure. So in the thermodynamic 
limit the steady state is spatially uniform except in the 
two extreme cases p = and p = 1. A glance at the re- 
sults exhibited in Fig.|fJ]for the local field and in Fig.[2]for 
the global fields reveals the great disparity of the effect 
of these fields on the steady-state properties of Axelrod's 
model. 
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FIG. 6: Normalized number of domains in the one- 
dimensional model for the local field as function of the chain 
size L. The filled symbols are for p = 1 (circles) and p — 
(triangles) whereas the unfilled symbols are for p = 0.99. The 
solid curve is 1/L and the other parameters are F = 5 and 
q = 10. 
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FIG. 7: Normalized number of domains in the two- 
dimensional model for the local field as function of the linear 
size L of the square lattice. The filled symbols are for p = 1 
(circles) and p — (triangles) whereas the unfilled symbols 
are for p — 0.9. The solid curve is 1/L 2 and the other param- 
eters are F = 2 and q = 8. 



Figure [7] shows the results of the local field in the two- 
dimensional Axelrod model. These simulations are in- 
credibly time consuming for p 6 (0, 1) due to the ten- 
dency to homogenization of the agents' cultural traits 
and the constant need to update the local fields. Hence 
our simulations are restricted to L < 200 and 10 3 samples 
only. As in the one-dimensional case, the simulation data 
in this range of p is practically indistinguishable within 
the numerical error and so we present only the results 
for p — 0.9. However, the data do not fall on the curve 
1/L 2 (solid line in the figure), which would signal the 
existence of a single domain. Rather, the scaled number 
of domains seems to go to zero much slower than L~ 2 
as the lattice size increases. This means that the aver- 
age (non-normalized) number of domains grows with L a 
with a < 2, i.e., it is non-extensive in presence of the 
local field. Of course, the small lattice sizes as well as 
the reduced number of samples used in this analysis does 
not allow us to make quantitative claims on the scaling 
laws for large L. We note that it was this difficulty to 
make inferences on the effect of the local field (even for 
a more manageable parameter set than that used in Ref. 
[ill ]) that motivated our analysis of the one-dimensional 
model. Finally, we note that the absorbing configura- 
tions are more heterogeneous for p — 1 than for p = in 
contrast to our findings for the one dimensional model. 



VI. CONCLUSION 

In our effort to clarify the effects of spatially uniform 
fields (interpreted as mass media) as well as of nonuni- 
form fields on Axelrod's model we have unveiled several 
interesting and unsuspected features of this well studied 



model of social influence. 

First and foremost, we have debunked the claim that 
in the thermodynamic limit the field-free heterogeneous 
regime becomes homogeneous in the limit of vanishingly 
small fields, i.e., p — S> regardless of the type of the field 
11]. More pointedly, we find that in the presence of a 



local field the steady state does indeed become homoge- 
neous but for all p £ (0, 1). This is hardly a surprise since 
the local field implements a majority- vote rule among the 
nearest neighbors of the target agent and so, when allied 
to the homogenizing interaction rule of Axelrod's model, 
it constitutes an insuperable force towards homogeniza- 
tion. However, we find that the system remains hetero- 
geneous in the presence of the spatially uniform fields 
which were originally introduced in Axelrod's model to 
study the effect of mass media on opinion formation Q . 

Second, we find that the spatially uniform but time- 
varying consensus media introduced in [8[ and the static 
media yield the same results in the thermodynamic limit. 
It seems the reason they produce different outcomes for 
small lattices is the correlation between the consensus 
field and the agents in the random initial configuration. 
This correlation becomes less pronounced for large lat- 
tices since even after application of the consensus rule 
there will be a rough balance between the values of the 
traits of a same entry of the feature vector. From the sta- 
tistical mechanics perspective, this is an important result 
because it allows us to make inferences about a more real- 
istic media type using the much easier to simulate static 
media. 

Third, we find that the heterogeneous absorbing con- 
figurations at p — and p — 1 are unstable. This insta- 
bility is reflected by the discontinuities that take place at 
those values even for finite lattices. Otherwise, the mea- 
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sures used to describe the steady state - average number 
of domains and size of the largest domain - are contin- 
uous functions of p in the range p G (0,1). We note 
that the field-free homogeneous configurations for finite 
L are stable in both one and two dimensions. However, 
in the thermodynamic limit they become unstable in one 
dimension (lower panel of Fig. [T]) but remain stable in 
two dimensions Q. 

Fourth, we find that the one-dimensional version of Ax- 
elrod's model Q yields essentially the same results as the 
more popular two-dimensional version. In particular, our 
findings about the impact of the three field types on the 
field-free heterogeneous regime of the two-dimensional 
model were corroborated by the one-dimensional model. 
Although this model exhibits some peculiarities, which 
were properly highlighted in the text, it can serve as an 
exceptional guide to our understanding of features that 
are difficult to unveil in the two-dimensional version, such 
as the effect of local fields. 

The characterization of the absorbing configurations of 
Axelrod's model in the thermodynamic limit, regardless 
of whether or not in the presence of a field (media), re- 
mains a challenge to statistical mechanics. In fact, we 
do not know much about the location and the nature of 
the phase transition in the (q, F) space @,|J], due mainly 
to the huge relaxation times the dynamics needs to set- 
tle in a homogeneous configurations, which grows as N 2 
where N is the number of lattice sites . In that sense, 
the appearance of dubitable claims about the behavior of 
Axelrod's model in the thermodynamic limit should not 



be surprising. This situation is worsened by the pres- 
ence of a field since lattices of intermediate sizes exhibit 
somewhat perversely a regime distinct from the thermo- 
dynamic one (see Fig. [T|). 

Nevertheless, we think that we have reached by now a 
good qualitative understanding of the effects of spatially 
uniform (mass media) as well as nonuniform (local) fields 
on the statistical properties of the two-dimensional Ax- 
ehod model. In the field-free homogeneous regime case, 
the presence of a vanishingly small field (i.e., p <C 1) leads 
to the breaking of the single mono-cultural domain into 
a macroscopic domain (media) and a multitude of mi- 
croscopic domains which occupy a finite area of the lat- 
tice [9(. The scenario is different in the one-dimensional 
model where the field pulverizes the giant domain into 
microscopic domains as shown in the lower panel of Fig. 
[U Regarding the field-free heterogeneous regime case, 
which was our main concern in this paper, the presence 
of a spatially uniform field does not produce cultural ho- 
mogenization, regardless of the field strength p. In con- 
trast, a local field that implements the majority-vote rule 
among the neighbors of the target agent leads to a ho- 
mogeneous steady state for p < 1. 
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